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We consider the two-dimensional (2D) flow in a flat free-slip surface that bounds a three- 
dimensional (3D) volume in which the flow is turbulent. The equations of motion for the two- 
dimensional flow in the surface are neither compressible nor incompressible but strongly influenced 
by the 3D flow underneath the surface. The velocity correlation functions in the 2D surface and in 
the 3D volume scale with the same exponents. In the viscous subrange the amplitudes are the same, 
but in the inertial subrange the 2D one is reduced to 2/3 of the 3D amplitude. The surface flow 
is more strongly intermittent than the 3D volume flow. Geometric scaling theory is used to derive 
a relation between the scaling of the velocity fleld and the density fluctuations of a passive scalar 
advected on the surface. 

PACS numbers: 47.27.Gs, 47.53. -Hn, 92.10.Lq 



I. INTRODUCTION 

We consider flows in a flat two-dimensional (2D) sur- 
face that bounds a three-dimensional (3D) volume with 
turbulent fluid motion. The boundary conditions are that 
of a free-slip surface so that the normal velocity compo- 
nent vanishes but the parallel components are not further 
constrained. To some extent, this is the situation of sur- 
face currents on a river or the sea, if waves and ripples are 
absent or can be neglected. Particles floating on the sur- 
face reflect the properties of the flow and provide an easy 
visualization. These flows have an obvious connection to 
oceanographical applications but they apparently 

have not been studied in further detail. Even in the re- 
cent theoretical and experimental investigations of the 
statistical properties of the particle distribution by Ott 
and co-workers [^|^ the modelling was based on random 
dissipative maps and not on the underlying flow. Simi- 
larly, Saichev and co-workers based their investiga- 
tion of passive particle advection and cluster formation 
on Gaussian random velocity fields, white in time. Thus, 
one of our aims here is to analyze the properties of sur- 
face flows arising from Navier-Stokes dynamics and to 
connect them to the statistics of particles floating on the 
surface, along the lines of our previous work on passive 
scalars advected in two-dimensional turbulent flows . 

The flow in the surface is two-dimensional, but it has 
properties that are different from that of the usual two- 
dimensional incompressible Navier-Stokes turbulence. 
Obviously, the velocity field is not constrained by mass 
conservation in the surface: there can be up- and down- 
welling motions in the (incompressible) bulk which on 
the surface will appear as sources and sinks for the ve- 
locity fleld. Velocity and vorticity can be exchanged with 
the bulk flow underneath, so that in the inviscid limit 
without forcing neither kinetic energy nor enstrophy are 



conserved. Such effects of compressibility arise also in 
experiments in two-dimensional turbulence in soap films 
and were discussed recently ]lO|~p^. 

The experiments of Goldburg et al. |jl^ are close to a 
laboratory realization of the kinds of flows that are in- 
vestigated here. A vertically oscillating grid in a tank of 
water is used to produce turbulence. If the water surface 
is sufficiently far away from the grid it remains essen- 
tially fiat and the surface flow can be visualized with 
mushrooom spores. The measured statistical properties 
of the flow are close to the ones that we will derive here. 
This opens the way to further experimental studies of 
the statistical properties of the velocity fleld and of the 
particle dynamics in this interesting flow. 

Finally, we should like to point out that the flows are 
also of interest from a theoretical point of view, since they 
can be thought of as flows with a symmetry plane: let the 
surface be 2; = and consider the reflection symmetry 
that under z goes to —2 the z-component of the veloc- 
ity field changes sign. This is a symmetry of the Navier- 
Stokes equation, that is to say, if initial conditions and 
driving preserve this symmetry so does the time evolved 
flow. 

It is our aim here to derive the equations of motion for 
such a flow (section IIA), to discuss the correlation func- 
tion if the 3D flow is turbulent (section IIB), to present 
numerical results on the statistics of the velocity, vor- 
ticity, and divergence flelds and on the boundary layer 
thickness (section III) and to derive a relation between 
the fractal dimension and the velocity correlation func- 
tion for the advection of scalars within geometric scaling 
theory (section IV). Concluding remarks are given in sec- 
tion V. 



'Present address: P.O. Box 208284, Yale University, New Haven CT 06520-8284 
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II. THE TWO-DIMENSIONAL FLOW IN A 
FREE-SLIP SURFACE 

In order to arrive at the properties of such a flow, two 
approaches are possible: one rehes on an expHcit repre- 
sentation of the flow fleld with proper boundary condi- 
tions and the other seeks to derive the equations of mo- 
tion from the 3-d Navier-Stokes equation. They provide 
complementary information on the system. 



u{x, y, z,t) = ^ U2n{x, y, t)z^'^ , 

oc 

v{x, y, z, t) = ^ V2n{x, y, t)z^'' , 



n=0 
oo 



w{x, y, z, t) = ^ W2n+i{x, y, t): 



.2n+l 



n=0 



(1) 

(2) 
(3) 



and 



A. Flow with refiexional symmetry 

We begin with the equations of motion and the effects 
of symmetry. Let u, v, w, be the x-, y-, and z-components 
of the velocity field (Fig. 1) and let p be the pressure field. 
The surface flow can be realized as flow in a symmetry 
plane, e.g. the plance z = if the velocity fleld is in- 
variant under the symmetry {u,v,w) — > —w) when 
z — z. 
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FIG. 1. Flow geometry. The surface flow lives in the 
shaded surface above a turbulent bulk flow, x and y are 
the coordinates in the surface and z is the one normal 
to it. In the numerical simulation the flow is driven by a 
shear flow in the cc-direction with variations in z. 

This suggests to expand the velocity components in 
power series in z, with only odd powers for w, and only 
even ones for u, v and p, 



p{x,y,z,t) = ^p2n(a;,2/,i)z^" ■ 



(4) 



n=0 



Substitution into the Navier-Stokes equation and order- 
ing with respect to powers of z gives for the two main 
components of interest, uo{x,y,t) and vo{x,y,t), the 
equations 

dtuo + (uo • V)mo = -dxPQ + lyAuo + 2vu2 + fu , (5) 
dtVQ + (uq • V)wo = -dypo + i/Avq + 2vv2 + fv ■ (6) 

The hat superscripts on position vectors x and R, on the 
velocity fleld u, on the gradient V and on the Laplace 
operator A indicate that they are restricted to the com- 
ponents X and y that lie in the surface. The driving of 
the turbulence is modelled by a volume force with com- 
ponents /„, and fw] as usual we expect that the statis- 
tical properties of the flow depend only weakly (through 
intermittency) on the kind of driving as long as it is con- 
flned to large scales, v is the kinematic viscosity of the 
fluid. 

For the full 3D velocity field mass conservation V • u = 
connects the normal and tangential components, viz. 
dzW = —dxU — dyV or, on the levels of the components 
[cf. i)-®], 



W2«+l = -{dxU2n + 9yW2„)/(2n -I- 1) , 



(7) 



Similarly, the pressure has to be determined from the 3D 
relation Ap = —V • [(u • V)u]. With the power series 
expansion from above this becomes to leading order in z. 



Apo(a;,y) + 2p2{x,y) = 2 [{dxUQ){dyVQ) - {dyUo){dxVQ)] - 2 {O^uq + dyVo) 



(8) 



Equations ^ , (^) and (^) are the equations of motion 
for the surface flow. Note that besides the surface ve- 
locity field {uq, vq) and the surface pressure po there are 
additional contributions from higher order terms in the 
power series in z: the viscous driving terms iyu2 and i>V2 
from shear effects in the normal direction and a contribu- 
tion {A)~^p2 to the pressure, also resulting from pressure 
variations in the wall normal direction. From the point of 
view of the flow in the surface, these terms are externally 
given and can hence be included in the volume driving 
forces. Note, however, that now the driving is no longer 



confined to large scales, as assumed in the usual scaling 
analysis. With all unspecified terms absorbed into effec- 
tive volume forces fu and fv , the equations of motion for 
Uq and Vq become finally 

dtUQ + (uq ■ VjuQ ^ -dxPo + v^uq + /„ , (9) 
dtVQ + (uo • V)t;o = -dyPQ + vAvq + fv , (10) 

The equations are completed by (^) with P2 = for the 
pressure. 

These equations have unusual properties. For instance. 
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dotting with u and integrating over a 2D volume, the en- 
ergy is not conserved in the Eulerian Umit where viscosity 
and driving are absent. With the local energy density, 



E{x,y,t) = {ul + vl)/2, 



(11) 



and using eq. (0) for n = the global energy balance 
reads 

dt{E{x,y,t))s = -{wi{x,y,t){E{x,y,t) + po{x,y,t))) s , 

(12) 

where (■)s denotes the average over the surface S. Thus 
energy is permanently put in and taken out according 
to the gradients of the z-component of u and the pres- 
sure fluctuations. Over large time intervals one can ex- 
pect that a flow equilibrium with constant average energy 
is established and that the time average of the right hand 
side of (12) vanishes. It seems that the lack of energy con- 



servation on short times gives rise to larger fluctuations 
and larger intermittency corrections (see below and [p^ ) . 
A similar discussion applies to the vorticity, and will be 
given in section IIIB below. 



B. Direct representation of a stress-free surface 

The alternative approach mentioned above starts from 
an explicit representation of the 3D velocity field that 
takes the boundary conditions into account. Consider the 
Fourier expansion of the velocity field, 

u(x, y, z, t) — u-^ ^{t) exp(iK • x) cos(n7rz) , (13) 
vlx, y, z, t) — ^(t) exp(iK • x) cos(n7rz) , (14) 

K,ri 

y, z, i) = ^ ^K,n(*) exp(iK • x) sin(n7rz) , (15) 



w[x 



K, 



where the summation extends over all 2D wave vectors 
K = {K^, Ky) in the surface and all integers n. The sine 
and cosine terms take into account the stress-free bound- 
ary conditions at the top and bottom surface, 

dzU = dzV = ui = for z — Q and z — 1 , (16) 

Incompressibility requires 



= 0. 



(17) 



One advantage of this representation is that it quickly 
leads to a prediction for the two-point correlation func- 
tions. In the 3D case Kolmogorov scaling without inter- 
mittency gives for the inertial regime a decay of ampli- 
tudes |uj^^p oc \k'^ + (ri,7r)2|~"/3 0. In the surface, 
the 2D amplitudes are obtained by summation on n. This 
brings in a factor of K that compensates the one missing 



from the volume element, which is K dK in 2D rather 
than fc^ dk as in 3D. As a net result scaling of the cor- 
relation function does not change. However, the absence 
of the third component of the velocity field reduces the 
amplitude to two third of its three-dimensional value. For 
the second order structure function, defined as 



S2{R) = (|u(x + R)-u(x)p) 
we expect in the inertial regime 

S2{R)^\S2{R)^R^'\ 



(18) 



(19) 



where again the hat distinguishes the 2D surface from 
the 3D bulk structure function. 



III. NUMERICAL SIMULATIONS 

The numerical simulations are based on a nearly ho- 
mogeneous turbulent shear flow bounded by stress-free 
surfaces at z = and z = 1 as given in Eg. (fig) . The 
velocity field is decomposed as in Eqns. (|13D-(|15|) and 
the Navier-Stokes e quat ions are integrated using a pseu- 
dospectral method ||l^Jl^. The simulations were done 
for Taylor Reynolds numbers Re\ — 59, 79 and 99, cal- 
culated from the streamwise velocity component u, i.e. 
Re\ = y'1^g/[{{dxu)'^Y/'^v] with root mean square veloc- 



ity Ur 



) ' . The properties of the 3D bulk flow 



are included here only to the extent that they are needed 
for the comparison between bulk and surface; they are 
further analyzed in p^ . 

Kolmogorov length = (j/'^/e)^/*, velocity Ur, = 
{evY^'^, and time scales — {y/eY^"^ are calculated from 
the 3D energy dissipation rate in the surface, i.e. 



(20) 



where indices 1,2, and 3 correspond to x, y, and z, respec- 
tively. For Re\ = 99 this dissipation rate in the surface 
is about 40% of the value in the bulk. 



A. Structure functions of the velocity field 

Form factors in the middle of the cell and on the sur- 
face are determined from 114 statistically independent 
snapshots of the turbulent flow. We focus on the scaling 
of the n-th order longitudinal structure functions, defined 
as 



S^iR, zq) 



[u(i -I- R, zo) - u(i, zo)] ■ R/i?|") . (21) 



In the bulk and without intermittency corrections the 
second order structure function is expected to scale like 
i?2 in the viscous subrange and like R^^^ in the inertial 
subrange p^ . A comparison between bulk and surface 
structure functions is shown in Fig. 2 for Rex = 99. 
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FIG. 2. Second order structure functions S^{R,zq)Iv1^ 
for i?eA = 99 normalized with Vr^ = {evY^^ and e from 
(^). Data from the surfaces at = and zo = 1 
are indicated by diamonds and connected by continu- 
ous Unes. Data in the bulk were taken at zq = 1/2 and 
are indicated by triangles and dashed lines. The inset 
shows the local scaling exponents from an extended self- 
similarity (ESS) analysis. Estimates between the vertical 
dashed lines, where the exponents are reasonably con- 
stant, give mean scaling exponents of 0.69 in the bulk 
and of 0.71 in the surface. These values are indicated 
by horizontal dashed lines. The two structure functions 

coincide in the viscous subrange but differ in the iner- 
tial subrange. This difference is predominantly in the 
amplitude and not in the scaling exponents, and consis- 
tent with (|l^). A local scaling exponent can be defined 
as C(^) = dlogS'^(.R, zo)/dlog_R. Unfortunately, both 
structure functions do not show an algebraic scaling be- 
havior at intermediate scales between the viscous and the 
forcing scale range for the values of Rex achieved here. 
Therefore, we apply the extended self-similarity (ESS) 
analysis ||l^ to the data. A local ESS scaling exponent 
can be calculated by relating local scaling exponents of 
second and third order structure functions, 



dlog[S^{R,zo)] 



(22) 



The distance vector R is taken in planes of fixed zq. As 
shown in the inset in Fig. 2 the bulk data give a local 
scaling exponent of about 0.69, in agreement with other 
observations, but in the surface the local slope is larger, 
about 0.71. This difference is small but statistically sig- 
nificant. Local exponents, based on averages over planes 
parallel to the surface, show almost no variation in the 
center of the cell but a clear trend when approaching 



the surface. This is demonstrated in Fig. 3 for the devi- 
ations SDn^siR, Zq) = Dn^(R, Zq) — n/3 from classical 
K41-scaling for orders n = 2 to n = 6 for different zq. 
The plane zq = 1/2 defines the middle between both 
free surfaces. Data sets for two different Taylor-Reynolds 
numbers Re\ = 99 (a) and 79 (b) are shown. The transi- 
tion from bulk to surface behavior can be used to define 



a surface layer, as discussed further in section [IIP 
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FIG. 3. Deviations of the local ESS scaling expo- 
nent from the classical Kolmogorov scaling for different 
heights Zq of the averaging plane and for different orders 
n = 2 to n = 6. Part (a) for Re\ = 99, averaged over 228 
samples. The deviations Dn^siR, zq) ^re the mean ob- 
tained for scales R between I877 and 4I77, as indicated by 
vertical lines in the inset of Fig 2. Part (b) for Re\ — 79, 
averaged over 254 samples. Here the exponents are ob- 
tained for scales R between IO77 and 27??. The vertical 
lines in dicate the surface layer that is analyzed further in 
section 
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The increase in intermittency seems to be connected 
with an increase in fluctuations due to lack of incompress- 
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ibility and lack of energy conservation. It is in line with 
results for passive scalar transport in models with com- 
pressi ble G aussian random flows that are 5-correlatedJn 
and with direct numerical simulations 
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time 

of isotropic supersonic turbulence. Interestingly, in the 
latter case the authors also noted a strong difference to 
incompressible turbulence near the crossover to the vis- 
cous subrange. In their case vortex filaments of high in- 
tensity and narrow regions of strong negative divergence, 
due to small-scale supersonic shocks, appeared. In our 
situation it is the fluctuations due to normal shear and 
normal velocity components below the surface that have 
a strong effect near the crossover to the viscous subrange. 

In the viscous subrange the amplitudes of the structure 
functions agree, but in the inertial subrange the surface 
structure function is smaller by a factor of 2/3. In the 
previous section we explained the reduction in amplitude 
in the inertial range by the reduction in the number of 
active degrees of freedom or Fourier modes. In the vis- 
cous subrange this argument does not apply, since we 
absorbed many additional contributions to the equations 
of motion into the volume driving force. The amplitude is 
larger since these extra contributions also have to be dis- 
sipated, but it should not exceed that of a 3D structure 
function since they originally come from a 3D flow. So in 
the viscous subrange the reduction in dimensionality is 
not noticable and the structure functions coincide. 



B. Structure function of the vorticity field 



p2[ . Note that in incompressible stationary turbulence 
the second order velocity and the second order vortic- 
ity structure function are connected by an exact relation 
n{R) = 2e/i/ - AS2{R) (mi]. This holds true in two 
and three dimensions, and has additional terms if the 
flows are not incompressible. 
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FIG. 4. Vorticity structure function ^{R)/t^ ^ for the 
surface flow at Re\ = 99. The data base is the same as 
for Fig. 2. 



Another quantity of interest in 2D flows is the vorticity 
Ltj = dxV — dyU and the structure function, 

= (|w(x + R) -a)(x)|2}. (23) 

In 2D incompressible turbulence squared vorticity is an 
additional inviscid invariant and gives rise to an inverse 
cascade of energy. In 3D a vortex stretching term (a;-V)u 
is present that prevents a conservation of enstrophy. In 
2D and for the normal component of the vorticity this 
reduces to a normal gradient of the velocity field which 
by incomprcssibility is connected to the divergence of the 
flow field in the surface. Thus, for the 2D free surface flow 
the vorticity transport equation reads 

dtuj + {u-V)uj ^ -oj{V -uj + vKuj + . (24) 

Thus, the non-vanishing divergence of the surface flow 
provides a kind of additional vorticity forcing in 2D. Con- 
sequently, squared vorticity cannot be an inviscid invari- 
ant, and no inverse cascade develops. 

The vorticity structure function for the data under- 
lying Fig. 2 is shown in Fig. 4. It saturates for larger 
separations to a non- vanishing value. Non- vanishing vor- 
ticity fluctuations were also observed in the experiment 
p3[ and interpreted as an indication that the observed 
features are not connected with turbulent surface waves 
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FIG. 5. Probability density function (pdf) of the sur- 
face vorticity component for the flow field of Fig. 2. For 
comparison a Gaussian pdf fitted to the central part of 
the distribution is indicated as well (dashed line). 

The strong intermittency of the flow is also reflected 
in the probability density function. Fig. 5 shows that 
the probability density function of the vorticity deviates 
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from a Gaussian distribution and has the exponentially 
stretched tails that are typical for intermittent quantities. 



C. Divergence of the surface flow 

The property that distinguishes surface flows from in- 
compressible 2-d flows most clearly is the divergence of 
the flow, which does not vanish for the surface flow. Snap- 
shots of the flow field, such as in Fig. 6, clearly show the 
presence of sources and sinks. 



27T 



liililfMi 



, -r-7f£::}j«?S3i«.; f?ij7v;";:;:';i;? 



mMMmMMmmm 




FIG. 6. A turbulent velocity field in the free-slip sur- 
face flow for Re\ ~ 99. The upper panel shows a vec- 
tor plot of the components u and v in the surface at 
z/d = 0. The lower panel shows a vertical cut through 
the box at the horizontal line (y/d = tt) marked in the 
upper panel. Regions of rising fluid and sinking fluid in 
the lower panel can be connected to sources and sinks 
near the solid line in the upper panel. A vertical slice 

across the flow underneath the surface allows to connect 
them to up- and down- welling motions below the surface. 
The corresponding contour plot of the divergence of the 
surface flow (Fig. 7) shows randomly fluctuating patches 
of sources and sinks. In the mean the flow is divergence 
free, ((V • u)) = 0, but the root mean square value does 
not vanish. Formally one can define a compressibility fac- 
tor ill 



< C = Ail — J-i < 1 



(25) 



(|VuP> 

which relates the mean square divergence to the mean 



square velocity gradient. For the surface flow and using 
only the velocity components in the surface, this becomes 



C = 



(iVuP) 



(26) 
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FIG. 7. Divergence V • u of the surface flow in units 
of the Kolmogorov time r^. Data are the same as for 
Fig. 6. In our simulations for Re\ — 99 we find C « 1/2, 

in good agreement with the Pittsburgh experiments [p^ . 
The relation of the denominator in (|2^) to the energy 
dissipation rate (|2^) is given by 



(|Vup)5 + ((V.u)2>5 



z.(|Vu|^)5(l+C), 



(27) 




FIG. 8. Radially averaged correlation function of the 
divergence field for Re\ — 99. The decorrelation length 
Ldiv is indicated by the vertical dotted line. 

The mean extension of regions with similar divergence 
can be determined from the correlation function. 
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Cdiv {R) 



V • u(x) V • u(x + R) 



(28) 



i/Ai 



(30) 



This correlation function is shown in Fig. 8. The first 
zero of Cdiv{R) defines a decorrelation length scale Ldiv] 
in units of the Kolmogorov scale L^i^ « 25. This scale fits 
rather well with the size of the largest patches in Fig. 7. 
As a consequence, the term in Eq. ( |2^ ) that contains the 
divergence of the velocity field describes a driving force 
that can be expected to be confined to the smaller scales 
in the flow. 




FIG. 9. Fluctuations of the normal velocity compo- 
nent as a function of the position between the sur- 
faces. The fluctuations are normalized in units of the 
square root of the mean turbulent kinetic energy — 
{v?)v + + {w'^)v- A linear extrapolation from the 

surface up to the value in the middle gives a boundary 
layer thickness 6 = 0.2 for Re\ = 59 and 5 = 0.1 for 
Re\ = 99. This is about a factor of 10 larger than the 
values estimated from Eq. ( ^ 



D. Fluctuations of the vertical velocity component 

In section IIIA we already mentioned the variations 
of the statistical properties with distance from the sur- 
face. They allow to identify a surface layer in which the 
transition from bulk to surface behavior takes place. This 
layer is different from the ones near rigid walls and not 
connected to friction but rather to the suppression of ve- 
locity fluctuations in wall normal direction. Dimensional 
arguments allow to determine the layer thickness S from 
a balance between the turbulent transport of wall normal 
fluctuations into the boundary, 



In both cases the size of the velocity field is estimated by 
the root mean square average of the wall normal velocity 
fluctuations, Wrms- Equating the two expressions gives 



v/Wr 



(31) 



as an estimate of the thickness. This is compared with 
numerical data in Fig. 9 where the z profiles of the 
wall-normal fiuctuations Wrms for two values of Taylor 
Reynolds number are shown. The boundary layer be- 
comes smaller with increasing Re\, as expected. The ab- 
solute values for the thickness of the boundary layer can 
be read off from the data by linear extrapolations of the 
profile slopes at the surfaces. The intersection of these 
straight lines with the corresponding maximum values 
of Wrms were used to define the boundary layer thick- 
ness (see Fig. 9). This gives values for 5 that are about 
a factor of 10 larger than the dimensional estimate, but 
consistent with its scaling behavior. This estimate for the 
thickness of a surface layer agrees with the interval over 
which the scaling exponents in (Fig. 3) change from bulk 
to surface values. 



IV. PASSIVE SCALAR TRANSPORT IN THE 
FREE-SLIP SURFACE 



A. Time scales 



Experimental and numerical studies show that the par- 
ticles floating on the surface of a fluid cluster in regions 
with down- welling and avoid regions with up-welling mo- 
tion ||,§. The patterns that appear have huge density 
variations that are best described by fractal scaling ex- 
ponents. As an approximation to the dynamics of par- 
ticles we can study the advection of a scalar density on 
the surface: it differs from true particles in that it has 
no inertia (the importance of which can be reduced in 
experiments by sufflciently small and light particles, see 
however |25[|2^ ]) and that it can develop larger gradients. 

Allowing for the compressibility of the flow fleld, the 
equations for the scalar density thus are 



dt(t> + V ■ (u0) = DK(t> - 



(32) 



(u- V)u; 



(29) 



where D is the passive scalar diffusivity. The Prandtl 
number is Pr = v/D. The following discussion will be 
confined to the two-dimensional fiow, so that all gradi- 
ent, divergence and Laplace operators act on the two co- 
ordinates X and y only; the superscript will henceforth 
be omitted. Expanding the second term in the above ex- 
pression then gives the evolution equation for the scalar. 



and the viscous dissipation of such fiuctuations. 



dt4> = -(u • V)(t> - (V • n)(t> + DA4> + , 



(33) 
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The input in scalar density is needed in order to com- 
pensate the diffusive losses. Since the equation for the 
scalar is linear, the natural amplitude scale for (jj is set 
by its root mean square value, iprms = After di- 

viding by (l)rms all terms have dimensions of inverse time 
and the time scales involved can be used to character- 
ize the different processes. Several of these processes also 
depend on the length scale I over which they are studied 
and so we introduce length scale resolved characteristic 
times. All terms can be made dimensionless using the in- 
ner scales of the turbulent velocity field in the surface, as 
discussed above Eq. ([20|). Different estimates of e bring in 
factors of order 1 (see e.g. relation ( pTj ) and the remarks 
below Eq. (pO|). Again, we use the energy dissipation rate 

Vrf, and Trj. 

i) is character- 



which is given by Eq. (|2C|) to compose 
The advection term, (u • V)(/) in Eq. 
ized by the advective time scale 

Tadv{l) = l/u{l) = — - 




(34) 



Z/77)i/3 in Kol- 
31) contains the 



where the scale resolved velocity u{l) w 
mogorov theory. The next term in Eq. (| 
divergence of the 2D surface flow and acts like a source 
or sink for the scalar. Its time scale is denoted TdiviJ)- 
The compressibility factor relates the divergence of the 
flow field to the root mean square velocity gradient [see 
Eq. (p5|)] which is connected with the energy dissipation 
rate e and thus the Kolmogorov time r^. Numerical sim- 
ulations and experiment Q indicate r^ivil) ~ r^/C^/^, 
where C is the compressibility factor (^6|) with a value 
of about 1/2. The efficiency of diffusion clearly depends 
on so that the time scale Tdiff{l) for diffusive smear- 
ing is Tdiff{l) « P/D. Finally, we have the forcing time 
Tf — iprms/ fc/i, which again is independent of spatial res- 
olution. 

In any given range of length scales, the process with 
the shortest time scale can be expected to dominate. So 
starting from the smallest scales we expect for an incom- 
pressible fluid first a diffusion dominated regime, then an 
advection dominated one and finally the input dominated 
regime. A Batchelor regime for the scalar is found if the 
diffusive regime extends beyond the Komogorov length, 
i.e. Pr ^ 1. In an incompressible flow, C = and Tdiv 
is infinite, so that there is no influence from the diver- 
gence. In the surface flows studied here, the estimate for 
C indicates that Tdiv is very short, of the order of the Kol- 
mogorov time. This implies that the advective regime is 
suppressed and that the statistics of the divergence dom- 
inates. This, finally, explains why the properties of the 
hydrodynamic flow do not seem to matter too much in 
the analysis of the particle distribution on free surfaces 
and why Ott et al. could explain the experiments using 
random maps |^-||]. 



B. Application of geometric scaling theory 



In order to connect the scaling of the velocity fleld 
to the scaling properties of the scalar, we use geometric 
measure theory ||2^ and the scaling ideas developed by 
Constantin and Procaccia |28-30|. A further extension 
of their work allowed for a scale resolved and Prandtl 
number dependent analysis |31 ^. The basic idea of the 
approach is to relate the fractal dimension 6g of the pas- 
sive scalar concentration, i.e. the scaling exponent (with 
respect to R) of the Hausdorff volume H of the passive 



scalar graph G{B^^^) = {(x,(/))|x G B)^' , = (f>{x)} 

(2) 

taken over a two-dimensional ball Bj^ of radius R, to 
scaling properties of the underlying turbulent flow that 
mixes the scalar. The following discussion will focus on 
the new terms relevant to the current problem; more 
details can be found in the above mentioned references 
p8|-^ and our previous work 

The basic quantity to be calculated within geometric 
measure theory is the relative Hausdorff volume of a sur- 
face of the normalized scalar density </> — (j)/4>rms, as 
given by 



j(2) 



H{G{Bf)) 
V{Bf) 



R 




V-^Pd^x, 



(35) 



where V{B^^'') — irR^ is the volume of a two-dimensional 
ball with radius R. The scaling exponent of the first order 
scalar structure function and fractal dimensions can be 
related by inequalities, which for the analysis arc assumed 
to be sharp Using the relation (?!)A(/) = A(?!)2/2- 1 V(?!)p 
and the equation of motion (^3|), the gradient under the 
integral can be replaced by 



|V0P=^-i5(u.V)^^-l0^(V.u) + lA^^ + - 



With 



(u • V)02 = V' 



(u4 



(36) 



(37) 



the first term on the right hand side can be expressed as 
a sum of two divergences. When substituted under the 
integral in Eq. (BS) the Hausdorff volume becomes 



V{B^n^) 



< 



\ 



1 



2D 



[-V-(U02) 



D(j>r 



d^x. 



(38) 
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The four integrals are denoted Ii through I4 and an- 
alyzed separately. The analysis of the three integrals /i, 
I3, and I4 proceeds as in the previous applications to 
two-dimensional scalar advection In particular, ap- 
plication of Gauss' theorem and the Cauchy-Schwartz 
inequality connects the first integral to the longitudinal 
structure function of the surface velocity field S2 {R) , 



4> 



D 



Fa is the fiatness of the passive scalar, 



(39) 



(40) 



If the correlations of ip decay rapidly this is essentially 
the volume average For a Gaussian veloc- 

ity field, = 3. Experiments and numerical simulations 
indicate strong ramp and cliff structures in the scalar 
field and thus some deviation from the Gaussian distri- 
bution |]3^-|35||, implying a scale dependence of F^. How- 
ever, we here restrict ourselves to first and second order 
correlations where intermittency corrections to the clas- 
sical Kolmogorov-Obukhov-Corrsin scaling are expected 
to be small and work with a constant F^. 

Exploiting the statistical stationarity of the passive 
scalar dynamics the fourth term, which contains the driv- 
ing of the passive scalar, can be expressed as 



■d^x, 



i?2 



1 



d X , 



(41) 



where R ~ R/rj is the radius of the disk in units of the 
Kolmogorov length. Using Gauss' theorem and Green's 



formula it can be shown |^ that the third term is sub- 
dominant compared to the fourth, 



(42) 



h < 2V/4 oc R 



and hence can be omitted in the following. 

Finally, we come to the new term, I2, which contains 
the divergence of the velocity field. Application of the 
Cauchy-Schwartz inequality gives 



V • u) d^x , 




(V-U)2,2 



(2). 



(43) 



where the root mean square of the divergence is measured 
in units of the Kolmogorov time t,,. The derivatives in 
the divergence term can be estimated from above using 
((V • u)2) < (|Vu|2) [see {^]. In the following calcula- 
tion this bound is not needed and the divergence fiuctua- 
tions in Eq. (43) can be taken directly from the numerical 
simulations. 

Combining (^ , (^ , and (|4|) we arrive at an inequal- 

(2) 

ity for the fractal dimension Sg of the passive scalar 

~ 1 129(1 at a fractal dimen- 



(1) 



)r(2) 



graph and thus via S, 
sion for the constant level sets (po = (/'(x), 



1 < 



d Ini? 



^lnh{R), 



(44) 



with 



h{R) 




u)2)i/2 + Pr^i?2 



(45) 



S2 = S2 /v'j^ is the longitudinal second order structure 
function in units of the Kolmogorov velocity. Obviously, 
if the last two terms under the integral dominate, then 
h{R) w it and S^g'^ — 2, implying a surface filling dis- 
tribution of the scalar. For sufficiently large scales R or 
large Prandtl numbers this is always the case. On the 
other hand, if the second term with the velocity structure 
hmction dominates, say S2 ~ , then 5^g^ = 3/2-1-7/4. 
Thus for the usual Kolmogorov scaling, S2 ^ R^^^ and 

4^) = 5/3. 

The Prandtl-number dependence of the fractal dimen- 
sion can be studied using as input the velocity correla- 
tions functions from our numerical simulations. Besides 
Pr also the prefactor T^^jTj, a measure of the strength of 



the scalar driving, is a free parameter in ( p5[ ) . The results 
for two different values of t^/tj are shown in Fig. 10. For 

many values of Pr a fractal dimension < 2 is ob- 
served. If the term t,,/t/ becomes large, either because 
of a small Tf (strong driving) or large (weak transport 
to smaller scales), the fractal dimension approaches that 
of a space filling fractal, i'^^^ w 2. 

In the experiments of Sommerer a fractal dimen- 
sion between 1.28 and 1.43 (denoted D2) was found. 
We find these values only in the transitional region, be- 
fore the inertial range is developed. The observations are 
consistent with Eq. (^) since it only provides an up- 
per bound and the observed values are smaller, indeed. 
Further comparisons between experiment and theory, us- 
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ing e.g. measured velocity correlation functions in (|45|), 
would be more than welcome. Some experiments are in 
preparation [ p^ . 




1 .0 . 

10 100 




10 1 00 

R/r, 



ergy nor enstrophy are conserved quantities in the Eule- 
rian, undriven limit. Moreover, in addition to large scale 
forces that maintain the 3D flow the surface flow is driven 
by small scale perturbations that come from transversal 
pressure variations and local gradients in normal velocity. 
As a result the scaling properties of the flow are essen- 
tially that of 3D turbulence, with an energy cascade in 
the inertial regime. The scaling exponents of the veloc- 
ity structure function are slightly larger than those of 
bulk 3D turbulent flows, indicating larger intermittency 
effects. We also observed a 2/3 difference in the amplitude 
of the structure functions between surface and bulk in the 
inertial regime. Several of the observed characteristics of 
the surface flow are in agreement with the measurements 
of the Pittsburgh group . 

We have also discussed the scaling properties of a scalar 
advected by the surface flow and have identified different 
scaling regimes. It seems that very often the dynamics 
induced by the divergence of the flow field is the fastest 
process, and that the advective properties of the flow are 
subdominant. This might explain why random mappings 
could succesfuUy be applied to the modelling of the par- 
ticle distributions but a more detailed comparison 
between theory and experiment is clearly needed. 

An important characteristic quantity of the surface 
flows is the compressibility factor C. The numerical sim- 
ulations and the Pittsburgh experiment using a verti- 
cally oscillating grid indicate both C « 1/2. With a sta- 
ble stratification of the fluid below the surface that re- 
duces vertical fluctuations it might be possible to achieve 
smaller values of C. This should open up the possibility 
to study the effects of compressibility over a larger range 
of C, both in connection with the intermittency contri- 
butions to the scaling exponents and with the scalar dy- 
namics in surface flows. 
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FIG. 10. Fractal dimension S^^'' for passively advected 
scalars for different values of the parameter a = Ti-i/Tf 
and the Prandtl number Pr. The underlying turbulent 
velocity field is the surface flow for Re\ — 99 as shown 
in Fig. 2 and = 3. 



V. SUMMARY 

The surface flows studied here are intermediate be- 
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